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Abstract. 

I study the time-evolution of a particle prepared in the ground state of an infinite 
well after the latter is suddenly expanded. It turns out that the probability density 
\^{x, t)p shows up quite a surprising behaviour: for definite times, plateaux appear for 
which |^(a;, t)p is constant on finite intervals for x. Elements of theoretical explanation 
are given by analyzing the singular component of the second derivative dxx'^ix, t). 
Analytical closed expressions are obtained for some specific times, which easily allow 
to show that, at these times, the density organizes itself into regular patterns provided 
the size of the box in large enough; more, above some critical time-dependent size, the 
density patterns are independent of the expansion parameter. It is seen how the density 
at these times simply results from a construction game with definite rules acting on 
the pieces of the initial density. 
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1. Introduction 

This sliort paper is devoted to some strange dynamical aspects around a problem which 
is often presented as the simplest one in quantum mechanics, namely the infinite one- 
dimensional well, although such a point of view can be seriously questioned (for instance, 
what about the Heisenberg equations of motion for the infinite well?). Indeed, when 
going beyond academic elementary questions, this problem is not simple, and even 
turns out to be somewhat tricky, all subtleties obviously originating from the infinite 
discontinuities of the potential, which generates an infinity of bound states with an 
energy En increasing without limit like the square of the quantum number n. This 
immediately entails that the propagator involves Gauss series (the Jacobi "i^s-function 
being one very special case [1]), which are known to possess quite uncommon features; 
as an example, Holschneider [2j shows that when the coefficients c„ of the series are 
oc the sum is a self-similar function in a precisely defined sense. Here, I only aim 
to give a brief account of intriguing results, together with a far from being complete 
theoretical explanation and proof. 

To be sure and strictly speaking, infinite discontinuities can be discarded on 
physical grounds, but they conveniently modelize a situation where the depth Vq of 
the well is much greater than all other relevant energies, and where the space variation 
of the potential occurs on a length scale / much smaller than all the others. Be it said in 
passing, for this reason, the classical limit of the infinite well is not a trivial point, due 
to the fact that one should first properly consider simultaneously the two limits I — > 
and Vq +oo, in order to check whether they commute or not and, if they do not, to 
choose the physically relevant limiting procedure for the considered case (for another 
example, see |3], §1.6). 

The problem at hand is the following. Given that the particle (mass m) is initially 
in an eigenstate of an infinite well, the well is instantaneously expanded to a larger 
size: what is the subsequent evolution of such a prepared initial state in the enlarged 
well? Some aspects of this problem have already been studied [1] - [7j; here I focus on 
results which are absent of these works and, up to my knowledge, seem unquoted in 
the literature. Obviously, any possible connection with an experiment would first of all 
require a proper analysis of various time scales, in order to be sure that the following 
theoretical framework is relevant to the experimental device. 

Let us now enter into the specific problem and precise the notations used 
throughout. Taking, for the non-expanded well, V{x) = when < x < a and infinite 
elsewhere, the normalized eigenfunctions are: 





(1) 



and vanish outside this interval; the eigenenergies are: 




(2) 
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n is a strictly positive integer, whereas Ti is the smallest time-period of any time- 
dependent state built as a linear combination of the ipnS. 

From now on, I assume that, the particle being in the ground state ipi{x) of the 
infinite well of width a, the latter is suddenly stretched at some time taken as t = 0, 
increasing in size from a to Xa with A > 1. Since the initial state ^(x, 0) = ipi{x) 
is not a stationary state of the dilated well, t) has a non-trivial time-dependence 
and, among other things, expectation values of the observables which do not commute 
with the Hamiltonian at t > show up actual time-dependence. I will focus on two of 
them, namely the probability density p{x, t) and the density probability current j{x, t) 
defined as usual: 

p(x, t) = |v|/(x, , j{x, t) = t)d^,^{x, t)] , (3) 

where 3 denotes the imaginary part, p and j are related by the local conservation 
equation dtp + dxj = 0. A few results concerning the averages of the position and the 
momentum of the particle will be briefly quoted at the end of the paper. 

On the other hand, the expectation value of the energy does not change since 
no work is done on the particle when the well is expanded; this obvious physical fact 
will be analytically checked in due time. As for the variance of the energy, it vanishes 
before the expansion, but turns out to be infinite once the latter has been performed 
(see section [5]). 



2. Wavefunction at t > 



The eigensolutions of the expanded well are simply obtained by making a ^ Xa in 
formulas ([1]) and ([2]), namely: 

^\,n{x) = sin^ = -^ipni^) , (4) 

V Xa Xa ^yx X 



Ex,n = . (5) 



for < a; < Xa, and: 

1 

A2 

Note that if is an irrational number, the two spectra i?„ and Ex^ „ have no coincidence 
at all. The dilatation of the well lowers each eigenenergy, and yields an increased energy 
density (in infinite space, the spectrum is continuous). 

The resulting state at time t > 0, \l/(x, t), can be developed on the complete 
eigenstates {V^A,n}n, and has an expansion of the form: 

+ 00 

^(x, t) = ^c„e^^^-"*V'A,n(x) . (6) 

n=l 

Note that, as thoroughly discussed by Styer[8] in connection with the classical limit, 
it immediately results that the motion is periodic, with the period T = A^Ti, since 
the expansion of \l/(x, t) only contains integer multiples of the circular frequency 
uj\ = X~'^uji] as obvious on physical grounds, enlarging the well increases the period T 
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of the motion: for an infinite expansion, the motion is not periodic since, among other 
things, the wavepacket would spread out ad infinitum. Also note that the wavefunction 
at time t is given by a Gauss series, i.e. a trigonometric series with time-oscillating 
factors of the form e'"^'^*, as contrasted to e'"'^* in a Fourier series. This yields quite 
rapid and irregular variations in time, all the more when the series coefficients decrease 
slowly with n, which is the case here (see ([8])). 

As for the coefficients c„, they are found by writing down the initial condition 
\E'(x, 0) = 'ipiix), and are thus equal to the scalar products (V'A.nl'^i), namely: 

2 . TTX . niTx , 

Cn = — 1= / sm — sm — — da; ; (7j 
aV A Jo a Aa 

note that the integral actually runs from to a, since ipi{x) vanishes for any x greater 
than a. A straightforward integration yields: 



2A^/^ sin 



rnr 



Cn = :r ^ 2 ' (8) 

TT A^ — 77,"^ 

SO that the wavefunction at time t > can be eventually written as: 

T / \ iA V — ^ sin -r- ; mvx ;„2, , i 

^ X, t = -\/- ^^^^^ ' , 9 

n=— oo 

for < a; < Aa, it being understood that \E'(x, t) vanishes outside the enlarged well. 

For any given time t, \E'(a;, t) is a continuous function of x and of t; this is recognized 
from the fact that the coefficients c„ behave like n"^ for large n, ensuring that the series 
in (ini) is uniformly convergent. Obviously, this is not true for the x- or t-derivatives of 
\E'(x, t) (remember that the potential has infinite discontinuities). 

By construction, each exponential function e„(a:, t) = e^*-^""^'^^*-' satisfies the 
Schrodinger equation ihdtCn = —^dxx^n, so that e*^dxxCn — Cndxx^n = 0: as it is the 
case for any stationary state in one dimension, the probability current constant in space, 
dxjst{x) = 0. This entails that performing a term-by-term derivation of the expansion 
([9]) to get the formal expression of the current j{x, t) related to ^(x, t) can only generate 
singular terms, arising from the difference between the derivative of a function, and the 
series of the derivatives; these singularities turn out to be Dirac functions, which means 
that, for a given time, j(x, t) is a piecewice constant function of x. Several examples of 
this will be given in due time. 

Note that making t = in the RHS of ([9]) leads to the function equal to 
\plja sin(7ra;/a) for < a; < a, and equal to zero for a < x < Aa, since \E'(a;, 0) = "^xix): 
in view of the sequel and considering the whole interval [0, Aa], this allows to say 
(trivially at this point) that the initial probability density shows up a plateau with 
a vanishing value for a <x < Aa. From this, one concludes that the following equality 
holds true for any a; G [0, Aa]: 

. ■« ~f-00 • Th'K 

- ^^e'^ = a-x sm— , 10 
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where 6{x) is the unit step function {6{x < 0) = and 6{x > 0) = 1), as well as all 
the other equalities obtained by a term-by-term derivation; all of them can be 2Xa- 
periodized in x if needed. The important point to realize is that the series in the LHS 
of (fTOj) is identically zero for any x such that a < x < Xa. It turns out unnecessary to 
define the step function for x = 0, since all the corresponding terms are multiplied by 
functions vanishing at this point. 

As we will see, one strange thing (among others) is that the probability density 
at time t also shows up plateaux (but not always with a vanishing value), in other finite 
intervals [xk, Xk+i] at given periodic times; this can be figured out as the recurrent 
ghosts of the initial flatness on [a, Xa]. 

Also note that if A — > 1 (no change of the well), all coefficients go to zero, except 
for ci which equals 1, as it must be. More generally, if A is a positive integer no, the 
indetermination for c^o is left by setting A = no + £, and by taking the limit e — 0. For 
A = no G N*, one thus obtains: 

c = hm ^ ^ ^^2±f^ = —— (11) 

Also note from that ^(x, T - t) = ^i^x, t), so that p{x, t) = p{x, T - t): 
at times t and T — t the two density distributions coincide, but since the two 
wavefunctions are complex conjugate, the two corresponding wavepackets have opposite 
group velocities; for the same reason the current satisfies j(x, T — t) = —j{x, t). Other 
symmetry properties can be found by inspection of the series (Q; for example, one easily 
sees that for t = T/4, \l/(x, T/4) = — \l/*(Aa — x, r/4), namely that at a quarter of the 
period (or at three-quarter), the density profile is even with regards to the middle of the 
dilated well. Other relations exist when both the abscissa and the time are changed, for 
instance one has 

^(x, t + T/2) = -^(Aa-x, t) (12) 

for any x and t. As we shall see, such symmetries play an important role, in particular 
to get closed convenient expressions for density and current at some remarkable times. 

Since the initial state ipii^) is normalized to unity, so is \l/(x, t) at any time; this 
can be checked by a direct summation of the series ^^=1 1'^^P (^^^ Appendix A). 



3. Probability density plateaux and hints for a theoretical explanation 

The surprise comes when plotting the probability density p{x, t) at different times. Some 
examples are given in figs. [T]- [HI which show that for very special times, the probability 
density assumes constant values in some definite intervals included in [0, Xa]. As already 
said, these plateaux can be figured out as the echoes of the flatness of \l/(x, 0) with a 
zero height in the range [a, Xa] for x. 

The theoretical explanation of the existence of the plateauxlies on arguments which 
could be more firmly grounded if mathematical rigor were required. The basic idea is to 
use the second derivative d^x^^ix, t) as an indicator, since its singularities determine the 
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2,0 n .'-.jn 



1= 1.5 




Figure 1. Probability density |^'(a;, i)P when the particle starts from the ground 
state of the undilated well; here, A — 1.5. Each curve is labelled by the time t, T being 
the period of the motion (see the text). 




Figure 3. Same as fig.[T]with A = 5.5. 
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abscissae where t) can have a cusp. Indeed, let us assume for definiteness that Xq is 
an abscissa to the left of which \E'(x, t) is increasing, and is constant on the right. This 
means that the first derivative dx"^ has a negative jump at xq, entailing that the second 
derivative contains an additive singular term oc S{x — Xq) with a negative weight, S{x) 
being the Dirac function (remember that if a function f{x) has a jump A/ at x = xq, its 
derivative is f'{x) + AfS{x — Xq) = f'{x) + -Dsing/, where /' is the ordinary derivative). 
At such a singular point, the derivative d^"^ has a jump, so that, generally speaking, the 
density |\E'p shows up a cusp. Due to the general properties of the Schrodinger equation, 
singularities are indeed to be expected in the second derivative in the presence of infinite 
discontinuities of the potential; they merely reflect, on a quantum mechanical level, the 
jumps of the velocity of a bouncing classical particle. These singularities, located at 
X = and x = a at t = actually move about in the interval [0, Xa] as time increases. 

The second derivative of \l/(x, t) is obtained by a term-by-term derivation of the 
expansion Q); by writing -23T2 = 1 H — it can be recast in the form: 



dx'^ a 



--^ij{x,t)+Di^^<i/ , (13) 



where: 



^^-.^ = -j;rAh E sm-e--e- (14) 



-00 



is the only singular part of The factor in the first term in the RHS of ( fT3ll 

is recognized as —^Ei and comes from the (ordinary) Laplacian operator in the 
(time-dependent) Schrodinger equation. Now, having in mind the well-known Fourier 
expansion of the Dirac comb: 

+00 +00 

^ e^— = ^(^-^) ' (15) 

n=—oo k=—oo 

it is realized that Dgjj^g\& embodies Dirac functions whenever the series in fll4p contains 
an infinite countable set of terms of the kind e^^™*''serx27r^ each having a coefficient which 
is independent of the dummy summation label. In order to explore this possibility, I 
rewrite the expression 0141) as follows: 

+00 

' ra=— 00 

First of all, note that for x = a, the first series in f fT6|) reduces to J2n=-oo^~^"'^'^^^^ 
generates a Dirac comb whenever e~'"^'^^* = 1 for an infinite countable set of values for 
n; this is the case if t = {p/q)T with p and q integers: for all values of n of the form 
kg {k integer), one has n^ujxt = k'^qp x 27r, which of the desired form: integer x 27r. At 
this stage, and considering only the first series in f|T6l) . it is seen that a cusp can occur 
for ^'(x, t), with {dx'^)a+ — {dx'^)a- > since the weight of 6{x — a) is then clearly a 
positive quantity. Note that the same argument also holds for all the points of the form 
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2^ = even integer but all the corresponding abscissas are outside the relevant interval 
[0, Xa] and may be ignored. 

This tells us that x = a is a good candidate, but this is just the beginning of the 
story, due to the existence of the second series in (fT6l) . To show what can happen, let 
us set X = a in both exponentials; the whole series then writes 



In fact, it can happen that for all countably set of "good" values of the integer n, the 
two exponentials cancel each other, annihilating the possibility for the point a; = a to 
be a cusp. For definiteness, and as an example, let us go back to fig.[Il and consider the 
curve t = T/2 where the density clearly shows up a "normal" maximum at x = a. For 
this case, one has p = 1 and g = 2 in the above notations, which entails that the good 
values for n are n = 2s (s integer); then, the only non- vanishing factors (1 — e^C^"^^/^)^ 
are for s = 1, 2 (3), but for s = 1 (3) and s = 2 (3), they have opposite signs, so that 
the two related Dirac combs indeed have opposite weights, and the singularity at x = a 
disappears. Thus, the point x = a is not always such a remarkable point. 

Let us now show that other values of the couple (x, t) can define the edges of the 
plateaux, without trying to give an exhaustive catalogue of all these possibilities, just 
aiming at giving a few sufficient conditions for that. 

The second exponential term in f[T6|) is equal to 1 for any n if 

(x + a)/{Xa) = r/s, r and s integers, and if n is an even multiple of s; the constraint 
< X < Aa entails 1 < r/s < 1 + 1/ A. This being realized, the conditions for the 
time-varying factor e~'"^'^^* are the same as above, namely t must be a rational fraction 
of the period T: t = {p/q)T. 

One example of such a case can be seen in fig.[Tl where A = 3/2. For t = T/A 
{uj\t = 271/4:), p = 1, g = 4 in the above notations. Close inspection reveals that x = a is 
indeed a cusp, as well as x = a/2 (take r = s = 1); there is numerical evidence, and this 
is analytically proved below, that these points are in fact the edges of a plateau. Note 
that the signs of the Dirac combs can be reversed; for instance with (x + a)/ (Aa) = r/s, 
if r is odd and n an odd multiple of s, e^^^'^+^^^'^^^Z**) = gi(2fc+i)r7r _ ^^^g examples, see 
the curves t = T/5 and T/10 in fig.[Il for which the density increases on the left and to 
the right of the plateau). 

Obviously, the existence of cusps is just a necessary condition for the occurrence 
of the plateaux. In order to analytically demonstrate their existence, one must generally 
prove that between two so identified given cusps, the density is indeed constant. 
This seems to be a rather intricate and difficult mathematical problem; in this short 
preliminary paper, I just intend to demonstrate this in a few specific cases, hoping to 
give a more complete and general proof in a future article. 




(17) 



n=— oo 
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4. Closed expressions for specific times 

It turns out that for some definite times tk, closed expressions of the wavefunction 
\E'(a;, tfc) can be written down. 1 will here consider only the three cases t = T/2^+^ 
with = 0, 1, 2, before showing the existence of quite another strange phenomenon, 
namely the fragmentation of the wavepacket and the existence of regular patterns when 
A is above a characteristic threshold Ac, depending on the specific time considered. The 
basic idea is to play with the time phase factors appearing in the expansion ([9]), and 
to express \E'(a;, t^) as a linear combination of the known initial wavefunction taken at 
different abscissae Xj. The generalisation for times of the form {p/q)T [p and q integers, 
p < q) seems quite feasible, although it promises to be somewhat cumbersome as long 
as a more elegant method is not available. 

A first observation is the following; at half of a period {t = T/2), a simple glance 
at the series IQ allows to establish the following equality: 

^(x, T/2) = -^(Aa -x,0) , (18) 

which is just the symmetry relation (fT2l) for t = 0; now, since \l'(x, 0) is known (this 
is 4'i{x), which identically vanishes between x = a and x = Xa), the equality f|T8|) just 
gives a closed simple expression for the wavefunction at this remarkable time. In the 
following, I show how such a method can be used for the other times defined above. 

4.I. The case t = T/4 

The case A = 3/2 (see fig.[I]) and the spectacular plateau occurring for t = T/4 draws 
attention on this peculiar time. To start with and to introduce the method, let us 
analyze the things in details, but for any A. The clue is simply to realize that for this 
peculiar time, the time- dependent exponential in ([9]) is equal to 1 if n is even, and to — i 
if n is odd. This allows to write \E'(x, T/4) in the form: 

^(x, T/4) = S2,o(x)-i52,i(x) , (19) 

where the two (real) sums 5*2, o and 5*2, i respectively correspond to even and odd values 
for the summation index n. Now that all the time-factors in the RHS of ( |T9i) are fixed, 
it is tempting to compare such an expansion with \E'(x, 0), which has a quite simple 
expression; observing that: 

^(x, 0) = 52,o(x) + 52,i(x) , (20) 

and: 

^(Aa-x, 0) = -52,o(a;) + S2,i(x) , (21) 

The two sums 5*2, ^ can now be expressed in terms of \I'(x, 0), which is known, thus 
readily obtaining the sum of the series ([9]) at this time: 

^(x, T/4) = -^[e-'^/*e{a - x) sin — - 

a 

e+W4^(o_Aa + x)sin^^^^^-^] , (22) 
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an equality which yields the closed simple expression of the density for any A: 



2 TTx \ ■ 2 7r(Aa — x) 



ap{x, T/4) = e{a - x) sin^ — + e{a - Xa + x) sin^ ' (^3) 

with still < X < Xa. Let me now choose A = 3/2; from (j23|) . it immediately results 
that: 

{sin2 2^ , < X < a/2 
1 , a/2 < X < a (24) 

cos2 ^ , a < X < 3a/2 

which proves the existence of the plateau between a/2 and a in this definite case. Note 
that this density is built in the following way: take the initial density, cut it into two 
pieces in the middle, translate the right part to the right of the distance a, draw a 
horizontal line between the two maxima, and divide the whole by a factor 2. We will 
recover such rules below, showing that the density at some other remarkable times can 
be built by playing with pieces of the initial density. 

The expression ( I23l) is true at t = T/4 for any A, and has two clearcut behaviours 
according to A < 2 or A > 2. In the first case, the two 6 functions are simultaneously 
non-zero in the interval [(A — l)a, a], so that: 

{sin^ ^ , < X < (A - l)a 

sin' ^ + sin' , (A - l)a < X < a (25) 

sin' , " a < X < A 

This shows that for 1 < A < 2 but A ^ 3/2, the function at T/4 has no plateau; in fact, 
for all such values of A, numerical plots show that the latter does exist, but for other 
times and not located between the two simple values a/2 and a. Clearly the relative 
simplicity of the A = 3/2 case is due to the fact that A is a "simple" rational number. 

Note that the x-derivative of the density is equal to the real part 3fJ(\l/*9x\l/); one 
easily checks from the expression ([22]) that, for A < 2 (and still t = T/4), 3fJ(\l/*(9i:\l/) 
never identically vanishes in a finite interval. Also note that if \l/(x, T/4) as given by 
fl22]) is a continuous function of x (as it must be), its x-derivative is not, although it is 
devoid of Dirac peaks due to the cancellation of \l/(x, t) at each jump of the derivative: 
the density p, a continuous function of x, can indeed shows up cusps. 

For A > 2, the two intervals [0, a] and [(A — l)a, Aa] do not overlap; then, expression 
( !23|) says that the wavefunction identically vanishes at t = T/4 for any x G [a, (A — l)a]. 
Examples of this are illustrated in figs. [2] and [3l it is seen that |\E'|' vanishes between a 
and 3a/2 for A = 2.5, between a and 9a/2 if A = 5.5. Thus, for t = T/4 and A > 2, 
the wavepacket splits in two distant parts: the particle is fully localized in two intervals 
separated by a finite one; in each interval, the profile is the clone of the initial one, 
just divided by 2. This fragmentation into identical curves will also occur for t = T/8: 
then, I will find four identical well-separated clusters, provided that A is greater than 4 
(see section l^73|) . each of them being one-quarter of the initial density |'?/'i(x)|' properly 
translated. 
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2,0 n 




0,0 0,5 1,0 1,5 



Figure 4. Probability density (upper smooth curves) and density current (lower 
piecewise constant curves) for three very close times near t = T /A. Note the extreme 
variability of the current in space (except at exactly a quarter of a period). The 
oscillations near the jumps arise from numerical troncations of the series and are 
obviously related to some kind of Gibbs phenomenon adapted to a Gauss series. 



As explained above, the current probability density j(x, t) is a constant piecewise 
function, generally having jumps when the first derivative of the wavefunction is 
discontinuous; otherwise stated, the jumps of j(a;, t) also occur whenever the singular 
part Dgjj^g\E' contains a Dirac comb. This turns out to happen in many points of the 
interval [0, Aa], as seen in fig. HI where all the functions have been numerically computed 
from the series ([9]). These plots also show that there is not necessary a direct relation 
between the jumps of the current and the edges of the plateaux, and reveals the irregular 
variation oid^^ , which is not always clearly visible on the plot of the density, all the more 
since a cusp can occur only if 'Si{^*dx^) 7^ 0. Now, starting from ([3]) with \l/(x, T/4) 
given by ( l22il . a straightforward calculation yields the piecewice constant expression: 

j{x, T/4) = -9{a - x)9[x - (A - l)a] sinvrA . (26) 

Again, the situation is quite different for 1 < A < 2 and for A > 2. In the first case, the 
current vanishes for < x < (A — l)a and for a < x < Aa; in the intermediate interval, it 
assumes the constant negative value sinvrA. Due to the conservation equation, the 
two points a; = (A — l)aand the only points where, at t = T/4, the time partial 

derivative dtp is non-zero. As contrasted, for A > 2, the current vanishes everywhere: 
not only at this time the wavepacket is split off in two fully disconnected parts, but 
the current between both regions is identically zero since there the wavefunction strictly 
vanishes. 

I mentioned above that, due to the central role of the Gauss series given in ([2]), it 
is expected that all quantities have a rather rapid and irregular variation in time. Such 
a fact is illustrated in fig. [5], where the probability density and current are plotted for a 
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fixed X as functions of time (remember tliat for t and T — t the densities are the same 
and the currents have reversed signs). At first glance, j{x, t) even looks like a singular 
function; remember that j is given by a double Gauss series. For A = 3/2, one has the 
symmetry j(f , t) = j{a, f - t). 




-3,0^ , , , , , -4,()^ , , , , , 

0,0 0,1 0,2 0,3 0,4 0,5 0,0 0,1 0,2 0,3 0,4 0,5 



Figure 5. Probability density p and current j as a function of time at a: = a/2 (left), 
middle of the well before the expansion, and x = Sa/A (right), middle of the well after 
the expansion. 



3,0n 



3,0n 




0,0 0,1 0,2 0,3 0,4 0,5 0,0 0,1 0,2 0,3 0,4 0,5 



Figure 6. Same as fig. [S] for A ~ 5/2. x — a/2 (left), middle of the well before the 
expansion, and x ~ 5a/ A (right), middle of the well after the expansion. 



Surprises in the suddenly- expanded infinite well 



13 



4.2. The case t = T/8 

I shall here follow the same arguments as before, the situation being a bit more complex. 
I first introduce four sums ^ {k = 0, 1, 2, 4) corresponding to the values n = ip — k 
of the dummy summation variable in the series (Q. Now, inspection of the time phase 
factors shows that one has: 

2^-1 

vl/(x, T/2^+1) = J2 e"'''^^'" S,.,,{x) , (27) 

k=0 

which has now the form of a Gauss sum. For N = 2, this gives: 

^(x, T/8) = 54,o(x) - 54,2(x) + e-'-/^[54,i(x) + S.^six)] . (28) 

I now follow the same idea as before, trying to choose definite abscissae Xi such that 
the space factor in ([9]) compensates in some way the dephasing due to the time factor. 
By trial and error, it is seen that the proper abscissae Xi which allow to express the 
various sums in terms of the initial wavefunction "^{xi, 0) are Aa/2 ± x and 3Xa/2 — x. 
First note that the sum S'4, i(x) + 5*4, 3 (x) is simply equal to the known quantity 5*2, i(x) 
already introduced in subsection 14. 1[ As for the difference S4^^o{x) — S4^2{x), I find the 
following: 

< X < ^ : S,,o{x) - = ^[^(y + X, 0) - - X, 0)] , (29) 

y < X < Aa : 54,o(a;)-54,2(x) = ^[^{-^ + x, 0)-v^(^-x, 0)] .(30) 

Great care must be exercized when writing the relations between the sums S2N^k and 
the values \&(xj, 0) due to the fact that the equality (fTOj) only holds for < x < Xa: 
outside this interval, the wave function vanishes, altough this is not the case for the 
sums since they are 2Aa-periodic functions. 

The above results eventually allow to write the following closed expression for 
\l/(x, T/8) valid for any A (setting ^ = x/a for simphcity): 

2a*(x, r/8) = - 0/<(0 + - ^)/>(0 + 

e-''^/^[^(l-Osin7re-^^(l-A + Osin7r(e-A)] , (31) 
where the two functions /< and /> are: 

/<(0=^(l-^-Osinvr(e+^) + ^^(l-^+Osinvr(e-^) , (32) 

A A 3A 3A 

/>(0=^(l+2-0sinvr(e--) + ^^(l-y + 0sinvr(e-y) . (33) 

Note that small times give more cusps that larger times; numerical runs confirm that 
the initial two cusps propagate through the interval [0, Xa] and multiply at the very 
beginning of the motion, before reducing in number when the time gets closer to half of 
a period. 
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In order to illustrate these results valid for any A, let me again take A = 3/2; then, 
the above formula give for 2^/a'^{x, T/8): 

< X < - : -(1 + i) sin — , (34) 
a a . . nx ttx 

- < X < - : — 1 sm cos — , (35) 

4 2 a a 

- <x <a : —1 sm (2 — i) cos — , (36) 

2 a a 

5a 7TX TTX 

0- ^ X < — : — sin (2 — i) cos — , (37) 

4 a a 

5a 3a , ^ Tcx , , 

- < a; < — : - 3-i cos — . 38 
4 ~ ~ 2 ^ ^ a ^ ^ 

This respectively gives the expressions for the dimensionless density ap{x, T/8) in the 
corresponding five intervals: 

1 2 . 1 1 2.1 

- sm TT^ , - , - + cos vr^ - - sm 2tt^ , 

115 

- + cos^ TT^ + - sin 27r,^ , - cos^ nC^ ; (39) 

note the plateau for a/4 < x < a/2, and the cusp at x = a, all features which 
are apparent in fig.[7] where is plotted the density ap{x, T/8) using the preceding 
formula, and the analytical expression (13T1) for the other A values. I checked that they 
give the same density as that obtained by a numerical calculation using directly the 
expansion ([9]). 



2,0 n 




Figure 7. Probability density p{x, T/8) calculated from the closed analytical 
expression ([3T|) . for three values of A. Note the coincidence of the three densities 
for < a; < a/4, and between a/4 and a/2 when A = 1.5 and 2.5. The fact that the 
density is constructed with pieces of |5'(a;, 0)p is clearly visible for the case A = 3. 
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Coming back to the general A case, the expressions ( 15T]) - show that \l/(a;, T/8) 
a priori shows up cusps at the following abscissa, which I precisely define for further 
reference: 

Xi = a , X2 = |A/2 — l\a ,X3 = Aa/2 , 

= 9(2 - A)(3A/2 - l)a + 9{\ - 2)(1 + A/2)a , 0:5 = (A - l)a . (40) 

Quite remarkably, they are equally spaced, being located at pa/ 4 {p = 1, 2, 3, 4, 5) for 
A = 3/2; for A > 2, where (A — l)a and Aa/2 merge, the cusp at (3Aa/2 — 1) gets out 
of the interval [0, Aa], but the cusp at (1 + A/2)a comes in so that there is still 5 cusps, 
which all remain in the latter interval for any A (see fig. [9]). I will come back to this in 
the following subsection. 



2,0 n 



0.123 T 




Figure 8. Probability density (upper smooth curves) and density current (lower 
piecewise constant curves) for three very close times near t = T/8). 



The current can also be easily computed; I find: 

irh 



j{x, T/8) 



2^2 ma2 



ci±(0 sm — + C3±(0 sm 



(41) 



where the functions Cr± depend on the considered interval; for x < Aa/2: 



ci4e) = ^(i-0M(i-^-e)+^(i-^+0]+^^(i-^+0^(i-A+0(42) 



and: 



C340 = ^(i- A + 0^(1-^-0 



(43) 



For X > Aa/2, one has: 



A 3A 

ci+(0 = ^(l+2-0[^(l-e) + ^(l-A + 0]-^(l-Y + 0^(l-A + 0(44) 



and: 



3A 

c3+(0 = ^(i-0^(i- Y + 



(45) 
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All this shows that j{x, T/8) is a piecewise constant function, as it must be. The density 
and the current are plotted in fig.[8]from the (truncated) series OH]) for three close times 
near T/8; note again the rapid variation of the current. For A > 4, the current vanishes 
everywhere. 

4-3. Fragmentation 

One sees in figs.E] and [3], which both correspond to A > 2, that for t = T/4, the 
wavepacket is split into two symmetric parts at the edges of the allowed interval for x. 
This is true for any A > 2, as a consequence of (l23l) : then, the two intervals [0, a] and 
[(A — l)a, Aa] do not overlap, so that the density is non-zero only for < x < a and 
(A — l)a < X < a; the two corresponding peaks are identical in shape, each equal to the 
initial density simply divided by 2. It thus turns out that for times T/4 (and 3T/4), 
the particle is fully localized into narrow domains and cannot be found between them. 
It can be said that, provided the expanded well has a size large enough, namely greater 
than 2a, there is the possibility for two identical bumps of width a localized at the edges 
of the box, with no density at all in between. 



The same phenomenon occurs for t = T/8 (and 7T/8): for A > 4, the density 
shows up four identical peaks, each of width a. Two of them are at the edges of the 
interval [0, Aa], the two others are on each side of the middle of the box. Interestingly 
enough, the onset of the four peaks occurs at A = 4, a threshold at which the cusps are 
equally spaced (two couples of them are degenerate because here A — l = l + A/2 and 




1 2 4 



Figure 9. Abscissas of tlie cusps as a function of A. The black segments show 
the domains where the density is non-zero; the hatched ones those where the density 
vanishes. Note that when A is above the threshold Ac = 4, the domains of non- vanishing 
density move away one from the other, but keep the same size and shape. 
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A/2 — 1 = 1). Once this has happened, the two middle peaks ("twin peaks") remain 
at the fixed distance a one from the other when A increases, being locahzed between 
Aa/2 ± a (central cusps), while the two edge peaks also remain unchanged and are still 
located between and a, and (A — l)a and Aa as A varies (see fig.[IO]). It thus turns out 
that for A above the critical value Ac = 4, the cusps dehneate the regions of vanishing 
and non-vanishing density: Aa/2 and (A/2 ± l)a for the central clusters, a and (A — l)a 
for the ones localized near the boundaries of the box. Again, one can say that when 
the size is large enough (now greater than 4a) , four identical peaks of width a can take 
place as indicated, and are independent of the expansion parameter A. 

To sum up this discussion, it can be stated that as far as A is greater than 4, 
the density p{x, T/8) is simply obtained by translating several times the initial density 
|'?/'i(x)p = p{x, 0) according to the formula: 

pix, T/8) = ^ 5^ P(x - L, 0) (A > 4) , (46) 

where the location of the maxima 1^ are a/2, (A ± l)a/2, and (A — l/2)a. Increasing 
the expansion factor does not alters the profile of each peak; the twin peaks stay locked 
around the center of the box, whereas the edge peaks are getting more and more far 
away. Remember that above this threshold, the current identically vanishes. 

1,0 n 

0,8 - 



0,6 - 




Figure 10. Probability density p at t = T/S, for A = 8 above the critical value Ac = 4. 
The fragmentation has occurred; the peaks now remain unchanged in size and shape 
when A varies and are located at the edges of the box, and on either side of the middle. 

Gathering the above the results with those obtained in the T/4 case, one can 
anticipate that for all times of the form T/2^+^, there exists a threshold Ac = 2^ 
above which a fragmentation into 2^ peaks occurs. The density profile consists of 
the elementary pattern jj^H'^^x, 0)p + I'i'i^^ - x, 0)p], and its 2^"^ - 1 clones 
translated by integer x Aa/2^~^; this is yet to be analytically proved in general, but 
numerical calculations allow to be convinced that this is true for any N (see fig. [TT] for an 
example). All this also confirms that many cusps exist at first times of the T-periodic 
motion, but remember that the time unit is precisely the period T = A^Ti, so that 
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tN = = A2Ti/2^+i > large does not mean small times: clearly, 

the two functions \['(a;, and \E'(a;, 0) have no resemblance, although the latter 

allows to build the former according the above rules. 

The above conjectures are done in the continuity of the analytical results given in 
this paper. Many other statements can be claimed in view of numerical evidence, but 
they still remain to be proved; let me give a few of them: 

(i) For all times of the form tu = T/M, M integer, there exists a threshold Ac(M) 
above which complete fragmentation occurs. 

If M is even, Ac(M) = M/2 and one gets a pattern of M/2 peaks located as above. 
If M is odd, fragmentation starts up at Ac = M, with M peaks; all peaks appear 
in twins except one, located near the origin. 

(ii) Fragmentation also takes place at times pT/M, with p integer. The number of 
peaks depends on whether p and M have common divisors or not. For instance, 
with M = 12, A > Ac = 6, one finds six peaks if p = 1, 5, three peaks if p = 2, 4, 
two peaks if p = 3, and a single peak at x = Aa if p = 6 (half-period). 

The method presented in this paper should be still efficient for proving these (and other) 
statements, although a more elegant procedure is highly wishable in order to make the 
analysis less cumbersome and more systematic. Work in this direction is in progress. 

0,3 n 



0,2- 




Figure 11. Probability density p bX t = T/2^+i avec = 4, for A = 20; here, the 
critical value is Ac = 2^ = 16. 

5. Other results 

After having focused on these rather outstanding behaviours, let mc take the oppor- 
tunity to add a few things for completeness, some of them being, as far as I know, 
unquoted in the literature. 
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As a first by-product, one can compute the probability P„(t) to find the energy 
Ex n when achieving a measurement of the energy at a time t > 0; according to one of 



the postulates of quantum mechanics, one has P„(t) 



|(^A,n|^(t))P = |c„P, namely: 



(A2-n2)2 



sm 



(47) 



if A is equal to an integer uq, the probability P„q can be found from ffTTj) . which yields 
Pno = I/uq. When A > 1, the distribution of the Pn is ever decreasing as a fonction of n; 
on the contrary, if A ^ 1, P„ has maximum for n ~ A, but the probability distribution 
is quite fiat (see fig. [T2l) . This maximum has a clear meaning on physical grounds: there 
is some kind of resonance in the vicinity of the two states having an energy Ex,n close 
to El, the initial (and constant) value for the average energy. It can be checked that the 
expectation value '^^gn* PnEx^n is indeed equal to Ei at any time (see Appendix A). 

p 

l,On n 



0,8 

0,6 



U= 1.5 




Figure 12. Probability distribution P„ for three values of A. 

Note that the variance of the energy is infinite, since the average {H'^) is given by 
a diverging series (P„ oc n~^, E\ „ oc n^). This is due to the fact that the prepared state 
effectively implies a large number of eigenstates V'A.n because the coefficients c„ have a 
slowly-decreasing algebraic n-dependence, so that high energies are relevant. This yields 
divergent energy fluctuations. 

The expectation values of the position, {x)(t), and of the momentum, {p)(t), also 
display interesting behaviour with time. An example is given in fig. [131 it is seen that 
the particle is periodically at rest on the average, since {x){t) is constant and equal to 
Aa/2 whereas {p){t) vanishes. This means that repeated measurements at those specific 
times would give exactly the same results as if the particle was in any stationary state of 
the dilated well. Measuring (independently) the energy would actually reveal the true 
nature of the state, giving for each measure one among all the possible energies Ex,n- It 
is also numerically observed that {x){t) is bounded by a/2 and (A — l/2)a: the particle, 
in the average, gets never closer than a/2 to the reflecting walls at a; = and x = Xa. 
The product AxAp is plotted as a function of time in fig.fT^ 

Note that the inverse process - sudden compression of the well, A < 1 - is 
impossible: one can not instantaneously generate a function vanishing for Aa < a; < a 
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Figure 13. Left: variations in time (0 < < < T/2) of the expectation value of the 
coordinate and of its variance. Right: same for the momentum. 

from a function which is finite in that interval. An infinite well can only be compressed 
with a finite rate; this case was analyzed in refs. [1] and [7] . 




0.5- X = 3I2 



' 0,0 0,1 0,2 0,3 0,4 0,5 

Figure 14. Variations in time {0 < t < T/2) of the product Ax Ap. 

As a final remark, let me mention that the limit A — > +oo can be achieved from 
the above formula, and indeed reproduces irreversible propagation in half-infinite space 
starting from the initial state ipii^)] the point is to observe that t) in ([H]) is a 

summation on the variable u = n/X, strictly equivalent to a Darboux sum, which quite 
naturally generates the Riemann integral over u in this limit (the differential element 
dz/ arises spontaneously from the factor 1/A in front of the summation). From (Q one 
can thus write: 

*(^.*)-!^/;;,.^«"-— -d.. (48) 

Note that the two zeroes of the denominator are just apparent singularities. Explicit 
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direct calculation allows to check that such an expression coincides with that obtained 
directly with the propagator of a free particle in R+i 

U{x, t; x', 0) = - / sin kx sin kx' e~'^^ dk {x, x' > 0) (49) 
Jo 

acting on the initial state ipilx) to build the state at time t according to the standard 
way, namely \l/(x, t) = U{x, t; x', 0)\l/(x', 0) dx'. In Appendix B, I show that the 
expression (148!) should not lead to misconceptions about the p-representation of this 
wavepacket. 



6. Concluding remarks 

As stated from the beginning, this paper just aimed to present a brief review of the 
rather strange results given above. Although the general existence of the plateaux 
is numerically established, I was able up to this point to give only some elements of 
theoretical explanation, and a genuine proof in the two particular cases t = T/4, T/8. 
Clearly, further investigation is required in order to provide a general demonstration, and 
also to define a systematic method for finding the precise points (x^, tfc) in space-time 
where such intriguing behaviour takes place. 

The fragmentation phenomenon also requires more attention; at this point, it can 
be conjectured that for t = T/2^+^, there exists a critical value Ac = 2^ above which 
spontaneous fragmentation occurs into 2^ peaks which are the translated replica of 
the initial density, divided by 1/2^; I gave an analytical proof only for = 1, 2, but 
numerical evidence allows to be convinced that this is a general result. For A = Ac, the 
density is an ordered finite lattice of adjacent bumps. It cannot be excluded that more 
complex patterns could be realized, going beyond the simple organization observed 
for t = T/2^~^^, although numerical calculations for times of the form pT/M {p and 
M integers) have, until now, unveiled spatial organization having the simple features 
described above. Last but not least, a transparent physical interpretation would be 
welcome, allowing to get physical insight explaining such amazing and counterintuitive 
behaviours. Work in these directions is in progress and, hopefully, will be published in 
the near future. 



Appendix A 

I here show how to check that the state "^{x, t) given by eq.Q is actually normalized to 
unity, and that the expectation value of the energy is indeed equal to Ei for any time, 
as it must be on physical grounds since no work is done on the particle when the well 
is suddenly expanded. 

Let us consider the function G{X, (p) defined as follows (A not an integer): 

2in0 
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this series is uniformly convergent for any real 0, so that G{(j)) is a continuous function. 
On the other hand, derivatives of G obviously contain generalized functions (the unit- 
step function and its derivatives). One readily sees that the definition ( l50i) allows to 
write: 

\mmm' = (|- o) - g(a, <p)]) . (si) 

\ / (/)=7r/A 

Let us now find G{X, 0), which is an even vr-periodic function of the variable 0. 
By differentiating twice the definition ( |50ll . one obtains a linear combination of the 
function G itself and of a Dirac comb. This means that the non-singular part of G 
precisely satisfies the differential equation S^^G + iX^G = for any g]0, 7r/2[; the 
general solution is Acos2A0 + i?sin2A0. The two constants A and B can be found by 
using the known equalities (Mittag-Laeffler expansions): 

+ 00 ^ 

n=—oci 

-\-oo / ^ \ ^ 

G(A, 7r/2) = V -^^^ = — ^ , (53) 

^ A2 - ^2 AsmvrA 

n=— oo 

which yield A = — ^ cot Avr and i? = ^, so that eventually: 

.co.A(2H-.) < < ./2) ; (54) 

A sm 71 A 

As anticipated above, G(A, 0) is a continuous function of 0, but its first derivative has 
a jump at = (vr), explaining the presence of the Dirac comb in the complete second- 
order differential equation for G{X, 0). Using now the rule expressed in (15T|) . one readily 
gets|(vi/(t)|vp(t))p = l. 

As for the average of the energy, one has: 

{H) = -^[G{X, 0)-G(A, vr/A)]- 

^Ir [G(A, 0) - G(A, 0)]^ ; (55) 



V^A 

the quantity in the brackets of the first line vanishes since it is proportional to 
"^{x = a, 0); due to fl5T]) . one is eventually left with 

{H) = E^\{^mm' = E, , (56) 

confirming that the expectation value of energy (H) is equal to Ei at all times negative 
or positive, as it must be. 

Appendix B 

I here intend to draw attention on a serious misconsception which could arise in view 
of the expression (HHj) . In order to make the discussion easier, I rewrite the latter as 
follows: 

^{x,t) = ^= e^W?i|,(p)e-ipV(2mft)dp , (57) 
V27r^ i-oo 
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where the function is: 

m = 'J^'^ipM , (58) 

where po = nh/a. At first sight, it looks obvious to claim that $(p) is the p- 
representation of the initial state, while the time-dependent exponential in the integral 
in (1571) is just the ordinary phase factor for the free particle starting in the $(p) state 
at initial time. As apparently trivial as it stands, this statement is simply wrong. In 
order to show this, let us draw a few consequences of it. 

First, it is easy to calculate the integral |$(p)pdp; one finds that it is equal 
to 2, instead of 1. Second, the true p-representation of the initial state can be easily and 
unambiguously calculated according to t = 0) = {27rh)~^^'^ e~'^^^^^ sin(7rx/a) dx, 
and turns out to be: 

1 3/2 

$(p,0) = -^^^(l + e-^-P/P°) ; (59) 

aside the fact that it comes out properly normalized to unity since \E'(x, 0) is, the 
function $(p, 0) is frankly different from the function $(p) given in fl58|) . Another 
drawback is that, due to standard rules of quantum mechanics for p-representation, the 
expectation value of the coordinate is: 

~ d ~ Yut ~ 

{x)(t) = ih dp$*(p)[— $(p) - -^$(p)] . (60) 
J-oo dp mh 

Since $(p) is an odd function of p, the integral vanishes, giving (x)(t) = 0, which 
is clearly incorrect: the wavepacket moves (and spreads out) in the free half-infinite 
space as time goes on. On the other hand, a non-vanishing integral would give a purely 
imaginaly expectation value since $(p) is a real- valued function, up to a constant phase. 

The error comes from the fact that everything stands in R+, instead of R. In 
other words, when a function f{x) arises as a Fourier integral of the form: 



hoo 



f{x) = I e'^''F{k)dk , (61) 

oo 



the equality holds true only for x > and one must not conclude at a glance (although 
this could happen to be correct) that the function F{k) is the Fourier transform of 
/(x): since all this holds true only if x > 0, and assuming that the Jordan's lemma is 
applicable, one can add to F{k) any function (f){k) which is analytic in the complex upper 
half-plane without changing the integral in the RHS of ( l6Ti) : the difference between $(p) 
and $(p, 0) is actually such a function (remember that ±po are apparent singularities). 
In other words, although the Fourier transformation /(x) F{k) is unambiguous, any 
intervening Fourier integral must be cautiously interpreted before to claim this is just 
the Fourier inversion formula; unconsidered intuitive identification can give incorrect 
results. Remind that for such functions defined in R+, the Laplace transformation is a 
much more secure method to proceed. 
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